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Abstract 

Dynamic deployment analyses of folded inflatable tubes 
are conducted to investigate modeling issues related to 
the deployment of solar sail booms. The analyses are 
necessary because ground tests include gravity effects 
and may poorly represent deployment in space. A 
control volume approach, available in the LS-DYNA 
nonlinear dynamic finite element code, and the ideal gas 
law are used to simulate the dynamic inflation 
deployment process. Three deployment issues are 
investigated for a tube packaged in a Z-fold 
configuration. The issues are the effect of the rate of 
inflation, the effect of residual air, and the effect of 
gravity. The results of the deployment analyses reveal 
that the time and amount of inflation gas required to 
achieve a full deployment are related to these issues. 

Introduction 

Ultra-Lightweight Inflatable (ULI) space structures have 
become attractive because they can meet structural 
requirements for space applications at a low cost. These 
thin membrane structures can be fabricated and 
deployed for millions of dollars less than conventional 
structures. The roles envisioned for inflatable structures 
include solar sails, space telescope sunshades, solar 
arrays, pressurized habitats in space, antenna reflectors, 
and optical telescope mirrors. Solar sails capture the 
momentum of sunlight photons. The area density of the 
available momentum is small. As a result, solar sails 
must be large. The thrust provided by these large sails is 
so small that the sails must also be thinner than paper to 
produce useful thrust vectors on the spacecraft. 
Billowing solar sails are less effective than flat solar 
sails due to the reduced momentum exchange provided 
by oblique incidence photons. To keep billowing within 
allowable limits, membrane tension is provided. The 
tension can be introduced by employing inflatable 
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spars to stretch the membranes, by spinning the sail to 
produce a centrifugal force, or by a combination of these 
methods. Three solar sail configurations [1,2] under 
investigation recently are shown in Figure 1. They 
include a square solar sail configuration, a disc solar sail 
configuration, and a heliogyro solar sail configuration. 
The square solar sail uses inflatable spars for 
deployment and for maintaining the tension in the sail 
membranes, while the disc and heliogyro solar sails spin 
after deployment to generate centrifugal tension forces 
which help maintain the desired configuration. The disc 
solar sail also has inflatable spars to assist its 
deployment. 

Developing and testing prototypes in space is 
prohibitively expensive. Deployment tests of solar sails 
in a laboratory do not accurately simulate their 
deployment in space. Even if a ground test is conducted 
in a vacuum chamber, the gravity effect can not be 
avoided. Gravity will affect the deployment dynamics, 
the structure's shape, and the regions experiencing self- 
contact during the inflation deployment. Also, 
deployments that involve spinning the entire structure 
are very difficult to conduct in ground based laboratory 
tests. Thus, laboratory testing of ULI structures should 
be supplemented with computational simulations. 

Many researchers have conducted computational 
simulations of the inflation deployment process for 
membrane structures. A closed form approximate 
analysis of the inflation deployment of a rolled (coiled) 
tube was derived by Steele and Fay [3], The model 
gives insight into understanding the unrolling process. 
However, the unrolling tube is supported by an infinite 
plane. The constraint of the plane limits the application 
of this model for simulating deployment in-space. Haug 
et al [4] and Salama et al [5] employed control volume 
(CV) and finite element methods to investigate folded 
space rigidizable antenna structures and folded inflatable 
cylindrical tubes, respectively. The CV approach, 
implemented in the PAM-CRASH [6] and LS-DYNA 
[7] codes, is based on an airbag inflation model 
developed by Wang and Nefske [8]. As this approach 
neglects the inertia of the inflation gas, it may not be 
adequate for simulating deployment with a high velocity 
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gas. To include the inflation gas inertia effect in the 
deployment simulation, the Arbitrary Lagrangian and 
Eulerian method (ALE) [9] needs to be used. However, 
the ALE method is computationally expensive since it 
requires the use of many Eulerian solid elements (for 
modeling the interior of the folded membrane structure.) 

Dynamic deployment simulations often employ explicit 
integration algorithms that require small time steps. 
Depending on the complexity of the structure and on the 
total time required for the deployment, on the order of a 
million time steps may be needed to conduct an inflation 
deployment simulation. It is known that increasing the 
mass flow rate of the inflating gas shortens the 
deployment analysis time. However, it is still an open 
question as to whether the fast inflation analysis results 
can be extrapolated to predict the structural deployment 
behaviors at slow inflation rates. 

Residual air, left inside the structure during packaging, 
can cause components to inflate simultaneously during 
the initial stages of deployment. In the NASA Inflatable 
Antenna Experiment (IAE) [10], unexpected 
deployment behavior was attributed to residual air. 
Thus, deployment simulations should include the effects 
of residual air. 

The purpose of this study is to employ the CV approach, 
available in LS-DYNA, to simulate the dynamic 
deployment of two inflatable membrane structure 
configurations, and to investigate issues that are critical 
to the success of in-space inflation deployments. The 
issues investigated are the effects of inflation speed 
(mass flow rate), residual air, and gravity. 

Deployment Simulations 
Inflation Modeling 

An inflatable membrane structure is modeled by 
discretizing it into a set of CVs. A longitudinal section 
of an inflatable tube and the inflation gas inside it form a 
CV. The boundary of the CV, called the control surface, 
changes its shape as the internal pressure is increased. 
Very thin and soft membrane elements are placed 
between adjacent CVs. These membrane elements 
stretch as the CVs deform. The total area of the 
membrane elements are referred to as the orifice area 
which controls the flow of inflation gas between 
adjacent CVs. A brief outline of this method is 
presented below. A more detailed description of the 
inflation modeling including contact algorithm is 
available in the literature, [7,8]. 


The incremental volume change for a CV depends on 
the net inflow-mass rate, the equation of state for the 
gas, and the dynamics of the membrane structure 
bounding the CV. Assuming all variables are known at 
time t-At, an approximation of the internal energy, 
E(t) , in the CV at time t is given by 

E(t) = E(t-At)+Cpm(t)AtTj„ (1) 

where is c p is the specific heat at constant pressure, 
At is the time step, T jn is the inflation gas temperature, 
and lh(t ) is the mass flow rate of the inflation gas. 

The gas mass density, p(t) , within the CV is 
approximated from the mass flow rate as follows, 

p(t) = [m(t - At) + m(t)At] IV(t- At) (2) 

where V(t-At)is the CVs volume at time t-At. 
According to the equation of state for an ideal gas, the 
pressure, p(t) , is calculated as, 

p{t) = {k-\)p(t)^- (3) 

m(t) 

where A' is the ratio of the specific heat at constant 
pressure to the specific heat at constant volume. 

The pressure is used as input to the finite element 
analysis to determine the structural configuration at time 
t . The equation of motion of the inflatable structure has 
the form 

[m]{d}+ [c]{d}+ [a]{Z)} = \R at } (4) 

where [m\, [c],and [a] are the global mass, damping, 
and stiffness matrices computed with respect to the 

current configuration, j R ext J is the external load vector; 
and {/)}, {d}, and {/)} are displacement, velocity, and 
acceleration vectors with respect to the current 
configuration at time t . Employing an explicit 
approach, the finite difference form of Equation (4) is 
expressed as 
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where At is the time step. Equation (5) is solved for 
{d} ( and the structure's shape at time t is then 
available. This method is very efficient when the 
damping is zero and the mass matrix is made to be 
diagonal by employing lumping procedures. Time steps 
on the order of 10" 6 seconds are typically required for the 
deployment models. As a result, the "wall clock" time, 
for a desktop workstation, can be on the order of weeks 
to simulate the deployment of a structure which in real 
time inflates in a few minutes. 

The work performed due to the volume expansion 
reduces the internal energy, therefore a modified internal 
energy, E(t) * , can be obtained according to the internal 
energy evolution equation, 
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Material Properties 

Both the coiled and Z-folded models are polyethylene 
tubes with a diameter of 3.82 inches and a thickness of 
0.006 inches. The values of the tube's Young’s 
modulus, Poisson's ratio, and density used are 25,000 
psi, 0.25, and 0.033 lbm/in 3 , respectively. The inflation 
is air at 70° F with a molecular weight of 28.97 
lbm/lbmole. The gas flows in at one end of the tube 
with an inflation mass flow rate (dm/ dt) of 0.1 xt 
lbm/sec. 

Coiled Tube Model 


The finite element model for the coiled tube packaged 
configuration is displayed in Figure 2. It consists of 
seven CVs and is created by employing a simple 
Archimedean spiral equation, r = ad + r . The symbol 

a represents a constant, 6 is the sweep angle, and r . is 
the initial radius as shown in Figure 2. The nodes are 
created starting at 6 = 0 and r = r and continue up to 

a user defined length. The unwinding of the coiled tube 
is shown in Figure 3. The volume and pressure as a 
function of time for each CV are shown in Figures 4 and 
5, respectively. The tube is fully deployed at time 0.35 
sec. The pressure is 8 psi, which induces a hoop strain of 
10 %. 


Z-Folded Tube Model 


where V ( t ) is the volume at time t computed using the 
divergence theorem. 

Baseline Studies 

The inflation deployment of a straight tube from its 
packaged configuration to its inflated shape is studied 
below. Two packaged configurations are modeled. In 
the first model the tube is packaged in a coiled shape, 
and in the second model it is in a Z-fold shape. We 
define the tube as deployed when it appears to be frilly 
extended and is free of kinks. Both models employ frilly 
integrated Belytschko-Tsay 4-node membrane elements 
[11], A number of simulations were made to observe 
the effects of the tube's finite element discretization on 
the opening dynamics. The studies indicated that the 
meshes employed below are refined enough to capture 
the physical differences occurring during deployment 
due to the rate of inflation, residual air, and gravity. 
Deployment simulations are made to address the effects 
of the rate of inflation, residual air, and gravity. 


The finite element model for the Z-folded tube packaged 
configuration is displayed in Figure 6. Inflated shapes 
of the Z-folded tube in various stages of deployment are 
shown in Figure 7. The Z-folded model contains four 
CVs with three orifices that are placed at the three fold 
lines. The time required to frilly deploy the tube is 0.19 
sec. The volume and pressure as a function of time for 
each CV are shown in Figures 8 and 9, respectively. 
The nominal volume of a frilly opened central CV is 91 
in'. The pressure at time 0.19 sec is 2.4 psi and it 
induced a hoop strain of 3.0%. 

Investigation of Deployment Issues Using Z-Folded 
Tube Model 

Effect of Inflation Rate 

Inflation mass flow rates ( dm/dt ) of O.Olxt lbm/sec, 
O.lx/ lbm/sec, and l.Ox? lbm/sec representing a slow, 
medium, and fast inflation are employed. Note, the 
medium inflation rate case is the baseline problem 
described above (see Figures 6 to 9.) 
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The slow inflation deployment sequence is displayed in 
Figure 10. It is difficult to fully deploy the tube with a 
slow inflation rate. The tube reaches a shape at 0.4 sec 
in which the last folded section will not swing upward 
until 1.38 sec. The sudden swing of the last section to 
obtain the inflated shape results in an abrupt volume 
increase of the last CV, see Figure 11. At that moment 
the pressure is about 10 psi (see Figure 12), which is 
much higher than the medium inflation rate case. The 
simulation indicates that the slow inflation process does 
not induce enough kinetic energy to flip the last section 
upward at an early time. This is similar to a folded 
garden hose that can not straighten itself by the internal 
water pressure. 

The results of the fast inflation rate simulation are 
displayed in Figures 13 to 15. As expected the fast 
inflation rate results in a quick deployment, full 
deployment achieved at 0.085 sec. The medium 
inflation rate requires 0.19 sec to deploy. Each section 
inflates sequentially for the fast inflation case as shown 
in Figures 13 and 14. The pressure buildups in the first 
and second CVs are higher than the buildups in the other 
two control volumes, see Figure 15. 

The orifices located at the connections between the 
control volumes control the amount of inflation gas 
going into each control volume per unit time. Two 
additional Z-fold deployment models are created by 
putting the orifices at different locations. In the first 
model, a five CV model is created by removing the 
orifices at the fold lines and installing orifices at the 
midsection of each of the four legs of the Z-folded 
model. In the second model, eight CVs are created by 
installing orifices at the three fold lines and at the 
midsections of each of the four legs. Deployments were 
simulated for each model using the three inflation rates 
mentioned above. The computed dynamic deployment 
characteristics were similar to the four CV Z-folded 
models discussed above. 

The required inflation time, the total mass of the 
inflation gas, and the pressures in the first and last CVs 
at the time of full deployment are listed in Table 1. In 
general, the deployment time is proportional to the 
inflation rate. The internal pressures at the moment of 
full deployment are nearly uniform for the slow and 
medium inflation rate cases (as indicated by the small 
pressure difference between the first CV and the last 
CV.) In the case of a fast inflation rate, the pressure 
difference between the first CV and the last control 
volume is large. Note, deployment with a low inflation 
rate requires a high inflation pressure and a large amount 
of inflation gas (to flip the last section of the Z-folded 
tube upward.) A fast inflation rate deployment also 


requires more inflation mass than the medium inflation 
rate deployment. It appears that the pressure distribution 
and the amount of inflation gas in the tube at the 
moment of full deployment can not be simply related to 
the inflation gas flow rate. 

Effect of Residual Air 

The effect of the presence of residual air on the 
deployment of Z-folded tube is investigated. Each CV, 
after the first one, is assumed to have a residual air mass 
of 2.703x1 (C lbm, which induces a pressure of 0.1 psi 
for a fully opened CV. The first CV is inflated by the 
inlet gas only and does not have any residual air present. 
At the start of the deployment simulation, the residual 
air expands each subsequent CV causing them to push 
each other apart. These initial motions allow the 
inflation gas to flow between CVs easier, and allow for a 
low pressure deployment. 

The results of the residual air deployment simulation at 
the slow inflation rate, O.Olxt lbm/sec, are shown in 
Figures 16 to 18. Again, the residual air causes the 
control volumes to expand simultaneously at an early 
time and push each other upward, see Figures 16 and 17. 
The residual air changed the deployment process to 
prevent kinking between the last two CVs and the tube 
reached full deployment at 0.470 sec, which is much 
faster than the slow inflation rate case studied above. 
The pressure of 1.5 psi at full deployment, see Figure 
1 8, is also much lower than the slow inflation rate case 
(pressure about 10 psi shown in Figure 12) obtained 
without residual air. 

Effect of Gravity 

Gravity adds a body force to the equation of motion. 
This additional force is expected to affect the 
accelerations, deformations and the contact between 
adjacent control volumes during inflation. Thus, the 
deployment process is expected to change. The gravity 
effect needs to be investigated before results from 
ground tests can be used to help predict deployment in 
space. 

The gravity effect is studied using the Z-fold model 
described above. The model shown in Figure 6 is 
modified by adding two 1.5 lb masses to the top end 
corner nodes of the tube (Nodes A and B). The masses 
are added to simulate an inertial resistance and gravity 
load similar to that expected when a complex structure is 
being deployed. The boundary conditions shown in 
Figure 6 are used again. A measure of the effect of 
gravity is the difference between the time to full 
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deployment with gravity enforced and the time for the 
case without gravity enforced. 

The results for the simulation without gravity enforced 
are displayed in Figures 19 to 21, and for the case with 
gravity enforced in Figures 22 to 24. Comparing 
Figures 19 and 22 a longer deployment time is required 
for the deployment with gravity, 0.60 sec versus 0.485 
sec. The intermediate deployment shapes are also quite 
different (see Figures 19 and 22). These differences can 
be observed by reviewing the expansion history of each 
CV. Comparing Figures 20 and 23, the third CV inflates 
faster in the presence of gravity. Also, a higher 
deployment pressure is required to overcome the 
additional forces due to gravity (compare Figures 21 and 
24). 

Concluding Remarks 

The control volume (CV) approach, as implemented in 
LS-DYNA code, was employed to simulate the inflation 
deployment of a tube when packaged first in a coiled 
tube configuration and second, when packaged in a Z- 
fold configuration. The CV approach is attractive 
because it uses a simple ideal gas law to compute the 
pressure change in each CV, and then uses the pressure 
to drive an incremental finite element analysis of the 
opening structure. Additional simulations were made 
with the Z-fold configurations to evaluate the effects of 
the rate of inflation, the presence of residual air in the 
packaged configuration, and the effect of deploying in 
the presence of gravity. 

The deployment simulations with different inflation 
rates indicate that a low inflation rate may not yield a 
desirable deployment sequence. With the low inflation 
rate, the control volumes may require unacceptably high 
pressures to open the fold lines (kinks) and obtain a fully 
deployed structure. It was also determined that the 
pressure and amount of inflation gas required for hill 
deployment are not simply related to the inflation rate. 
Deployment results for a fast inflation rate can not be 
scaled to predict the structural behaviors at a slow 
inflation rate. 

A very small amount of residual air can dynamically 
open the control volumes, pushing them apart from each 
other. This helps prevent the "kink formation" and 
results in a much shorter deployment time with a lower 
deployment pressure. 

Gravitational forces can impair the opening of the fold 
lines in a packaged tube. This is similar to the case of a 
deployment at a low inflation rate. Fligh deployment 


pressures (over those required in space) may be 
required. Thus, deployment test results obtained on the 
ground will not fully predict a deployment in space. 
Flowever, ground tests can be used to validate an 
inflation simulation model by correlating the 
computational results with the ground test data. The on 
ground validated model can then be used to better 
understand deployment in space. 
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Table 1. Inflation time, internal pressure, and total mass of the inflation gas for Z-fold deployment simulations. 


Inflation rate, 
dm / dt 
(lbm/sec) 

Inflation time *’** , 
r 

(sec) 

Pressure**, p( t) 
(psi) 

Total Inflation 
Mass**, 
m( t) 

(lbm) xlO' 3 

First 

CV 

Last 

CV 

Ave. of 
all CVs 

O.Olxt 

4-CVs 

1.38 

9.26 

9.23 

9.24 

9.52 

5-CVs 

0.98 

5.31 

5.27 

5.28 

4.80 

8-CVs 

0.75 

3.29 

3.25 

3.26 

2.81 

O.lxt 

4-CVs 

0.19 

2.30 

2.13 

2.24 

1.81 

5-CVs 

0.22 

3.11 

2.75 

2.86 

2.42 

8-CVs 

0.20 

2.57 

2.20 

2.33 

2.00 

1.0 x? 

4-CVs 

0.085 

5.27 

3.69 

4.50 

3.61 

5-CVs 

0.085 

6.47 

3.53 

4.41 

3.61 

8-CVs 

0.080 

6.09 

2.62 

3.75 

3.20 


* Time step, At =8. 35X10 -6 seconds 

** r = time at lull deployment 
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Figure 4. Volume of each control volume (CV) as a function of time. 
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Figure 5. Pressure of each control volume (CV) as a function of time. 


Figure 6. Finite element model of a Z-folded tube. 
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Figure 7. Z-foldec 
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Figure 9. Pressure of each control volume (C'V) as a function of time. 




t=0.10 sec 


Figure 10. Different deployment stages of the Z-folded tube under a slow inflation rate 
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Figure 13. Different deployment stages of the Z-folded tube under a fast inflation rate. 
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Figure 14. Volume of each control volume (CV) as a function of time for the fast inflation case. 
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Figure 15. Pressure of each control volume (CV) as a function of time for the fast inflation case. 


Figure 16. Different deployment stages of the Z-folded tube with residual air. 
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Figure 2 1 . Pressure of each control volume (CV) as a function of time for the no gravity case. 
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Figure 22. Different deployment stages of the Z- folded tube with end mass and gravity. 
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Figure 23. Volume of each control volume (CV) as a function of time for the gravity case. 
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Figure 24. Pressure of each control volume (CV) as a function of time for the gravity case. 



